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1 Summary of spectral measurements from li^tning. These data have beer 

normalized to a distance of SO km and unit bandwidth. The spectrum is 
proportional to (frequency)"* in the frequency range from a few lO’s of 
kilohertz to about 100 megahertz 21 

2 The current pulse used in the calculations and its spectrum (i.e. magnitude 

of its Fourier transform). The current pu'.se is representative of current in 
first return strokes. In this example, a » 26 ■ 4 x 10^ sec"* ; 0 * 8 x 10^ 
sec"* ; 7 ■ 10^ sec"* and Iq * 30 kA and I] * 2.5 kA 22 

3 Example of a tortuous channel produced by the computer simulation. . . 23 

4 The spectrum calculated from a straight channel (about S km long) and a 
tortuous channel (the one shown in Figure 3). Only the envelopes of the 
calculated spectra are plotted. The current pulse shown in Figure 2 was 
employed assuming a velocity of propagation equal to c/2. The observer 

is SO km from the channel 24 
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The calculated spectra (t'ruin Figure 4) with representative data points 
superiit'|H)sed. The data were taken from Kimpara (Kimpaia, 1965) 
normali/ed to 50 km and unit bandwidth 


THF INFLUENCE OF TORTUOSITY ON THF SPECTRUM OF 
RADIATION FROM LIGHTNING RETURN STROKES 

David M. Le Vine 
ABSTRACT 

An investigation has been made of the influence of tortuosity on the 
spectrum of radiation from lightning return strokes. The shape of the 
spectrum obtained by hicluding effects of tortuousity is in keeping 
with data: The spectrum has a peak in the correct frequency regime 
followed by an initial decrease as the inverse of frequency. This 
spectrum is in better agreement with data than the spectrum predicted 
by the same model without tortuousity (i.c. the long straight channel), 
which decays at a rate proportional to I 

These conclusions were derived from a piecewise-iinear “transmission 
line” model for the channel using in one case simplifying assumptions 
to arrive at closed form expressions for the spectrum and also from 
numerical calculations on simulated channels. Tlie analysis indicates 
an eventual transition to I/t'^ decrease even for the tortuous channel, 
which suggests a means for testing the model and gaining insight into 
return stroke current and channel parameters. 
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THE INFLUENCE OF TORTUOSITY ON THE SPECTRUM OF 
RF RADIATION FROM LIGHTNINC 


INTRODUCTION 

Knowledge of the spectrum of radiation from lightning is important both for under- 
standing the dynamics of the lightning flash as well as for assessing the RF interference 
environment during thunderstorms. Spectral measurements have been reported from fre- 
quencies of a few kilohertz (Horner, 1961; Taylor, 1963; Watt and Maxwell, I9S7) to 
frequencies near a GHz (Kosarev, et al., 1970; Hewitt, 1957) and several composite spec- 
tra have been deduced from these data by reducing the independent measurements to 
common units of distance and bandwidth (Horner and Bradley, 1964; Kimpara, 1965; Oh. 
1969; Cianos, Oetzel and Pierce, l972;Oetzcl and Pierce, 1969). These composite spectra 
have several features in common; It is generally agreed that the spectrum peaks near 10 
kHz and then decreases roughly inversely with frequency (i.e. as \jv) upto several hundred 
Mil/.. At higher frequencies the data is more scattered: Kimpara (Kimpara. I965)extn polates 
the Mv dependence to higher frequencies whereas Oh (Oh, 1969) infers from the data a 
change from Mv \o a decrease as approximately Xjv' in the decade between 100 MHz 
and I GHz. Representative data are shown in Figure I. 

Theory to rationalize these observations with contemporary best information on 
lightning parametei's have been less widely reported; and in fact, the measurements pose 
some problems for the modeller. For example, using as a return stroke model a long r traiglit 
channel (current filament) driven by a current pulse whose shape is consistent with measure- 
ments, such as the paired exponentials of the Bruce-Golde form or as modified to include 
continuing current (Bruce-Golde. 1941 ; Uman, 1969; Price and Pierce, 1974), and assuming 
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that the ciirrcnl piilu' propagates along the channel with velocity, v (i.e. l(/,t) ■ l(t-z/v)>, 
one has the “transmission line*’ model presently in use (Dennis and Pierce, I^M; Uman 
and McLain, l‘>f>9|. It is nut particularly dilTicult with this model to calculate the electric 
Held radiated to a distant observer, employing for example the I'raunhuffer approximation 
(e.g. Le V'ine and Meneghini, l*)7b). Assuming typical channel lengths and velocities of 
propagation, this calculation yields a s|H‘ctrum whose peak is in the correct frei|uency 
range (tens of kilohertz); however, the high frequency asymptote decays as l/»^, not in- 
versely with the first power of frt*t|uency as the data suggests. 

It is the purpose of this paper to show that by introducing tortuoi'sity into the channel 
model that is. using the transmission line model but allowing the channel to be irregular 
rather than straight — one obtains a s|>ectral decay proportional to I/p. llie spectrum 
obtained with tortuousity included also peaks near 10 kli/. and then decreases roughly as 
l/p until, at an upper frequency determined by the scale of tortuousity, it begins to de- 
crease as l/p2. The shape is similar to that reported by Oh (Oh, I ‘>69). 

Relatively little attention has been devoted to the effects of tortuousity on the fields 
radiated from lightning return strokes. Hill (Mill, I96P) has treated the subject using a 
scalar analysis based on the moment approximation and the Bruce-Golde form for motion 
along the channel. This restricts the analysis to VLF (tor which it was intended). On 
the other hand, at fret|uencies of 10 MHz or so, the free space wavelength of the radiated 
fields is comparable to the typical step size (the channel is formed in discrete pieces dur- 
ing the stepped leader processes; Uman, 1969) and the individual elements of which the 
channel is formed may become effective radiators. In this paper, the effects of tortuousity 
are analyzed by adopting a piecewise linear model for the channel and assuming that the 
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current pulke propagate* along the channel with conttant velocity on each tegment: 
T(r,t) ■ $ f(t>$>r/v). This U the “trantmittion line” model applied on a piecewiu* hakU 
to the channel. Both an analytical model for the tpectrum. based on limpliryiiig aktiimp* 
tiont, and numerical results will be presented. The numerical results are obtained by 
having a computer produce channel rcaii/ations in accordance with pre-assigned statir'ics 
for mean length and orientation of individual elements. Then the spectrum is computed 
for each realization using formulas for the radiation from each individual element keeping 
track of the phau* as the current pulse propagates up the ciiannel (Le Vine and Meneghini, 
197Kb). The analytical .esults are obtained by simplifying the expressions for radiation 
from each element sufficiently that sums and averages can be iHrrformed explicitly. The 


results agree. 
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MODM 

A 

('on»iiicr a filament L iinitk lung m the direction driven by the traveling current 
wave: 

Ttr.t) - ntt-«*r/v) 

where v it the velocity of propagation along the channel, and avtume that the filament it 
arbitrarily oriented alxive a conducting plane at / * U. 

One can obtain a volution for the fieldv privjuced by tuch a current filament by 
transforming to the frei|iieiK 7 domain and solving in terms of potentials. IX>ing to. one 
obtains the following form for the magnetic vector potential: 


Atfo. p) • n !(»') 


/, 


libmeni 


^jkrjV-T. j ^ 4- } jr, 

4»R 


rikK 

TJTr' 


(II 


where the supersc’ripl tilda (^) denotes a I'ouriei ransform (for example, f (ri is the 
l-ourier 'ransform of f(tl) and where R » distance from the source point T, 

to the observer fj, and R’ ■ distance from the image source point Tj ■ 

r, - 2/ir, • r) to (he observer, and (I ■ • /) and t) ■ c/v. Hie electromagnetic 

fields are obtained from I'liuation I by means of the formulas: 


I (F„, r'l » jkc ^AiFq, I') + V ( r • Alj 


ll(F„, I'l » — C * A 


In general, the integral in l-quation I can not be evaluated and approximations must be 
made. The I raunhoffer approximation (i.e. far field approximation! will bv* employed 
here, (In the s|vcial case v * c the integral can be evaluated exactly in closed form (Le 
Vine and Meneghini, |d7Sa; Schelkunoff. I d52): although for real lightning, v < c, and 
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generilly v is not even constant along the channel.) To obtain the FraunholTer approxi- 
mati jn. Equation I is substituted into Eo> <tion 2 and it is assumed that the filament length, 
L, is smaller than the distance, R, from filament to observer. Then a Taylor series expan- 
sion is made of R about the Hlament center at Letting • IFo~^cl assuming that 
L/Pf < I, kL^/pc « 2a, kp^ » 1 and keeping only lowest order terms in kp^, one ob- 
tains the following form for the radiated fields; 


IT(To,i») ■ y/fUi f(a) - (T'VPc^ ^®cl “ 

1^ - (^' * 1 * — , I' (•») I 

4* Pc ) 

^ A ci^Pc A ei^Pc 

f(p) ( VPc * ^ Z I(p)-(^Pc XX ) T“* • 

4»Pc 4irpj 


cjkT,(t-r,) ^ 


ejV4kL(fj-^-VPc) - e-j'“‘L(n-^VP, ) 


Tj -T* VPc 


jkL cjkn(^*fc) sine ('/ikL(r 7 -^*V. 0 £)l 

L »7 - X • V Pc J 

jkL sine ('/4kL(q4‘ • vPc')1 


Oa) 

(3b) 

(4a) 

(4b) 


All interesting interpretation of Equations 3 and 4 is obtained by noting that the 
phase terms (the arguments of the exponentials) can be written in the form j2apr where 
r »(Pc/c - $ *Fc/p) ± '/i(L/p - L/c VPc)- The first term in parenthesis is just the time 
required for a signal to propagate from the filament center to the observer plus an arbi- 
trary constant which is necessary if several elements are to be compared. The 
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remaining parcnlhctit is ll;c binomial correction of this time a>rre)»poiuling to propagation 

fiom the upper and lower ends of the filament it — phis the time required lor 

2c 

the pulse to propagate from one end of the filament to the other (L'Vl. Tliat is, within 
the limitations of the Kraunhoffer approximation the radiation may be regarded as consis- 
ting of two pulses which emanate from the ends of the filament but with a relative delay 
equal to the time required for the pulse to propagate the length of the lilament, Interes- 
tingly, this is exactly what one finds to be triK* for the exact solution when v ■ c (Le Vine 
and Meneghini, l<)7Ka). 

Ass4>cialing the radiation with the end points of the filament permits a simple interpre- 
tation of effects of tortuosity. When me current pulse begins to propagate up the filament, a 
pulse of radiation (having the shape of the current pulse) is radiated from the bottom of the 
filament, and nothing more hap|H*ns until the current pulse reaches the top of the filament; 
then an identical pulse, but of opposite sign, is emitted and the total radiation is now the 
sum of these two terms (Uman, Mclain & Krider, l‘)75: Le Vine & Meneghini, iy7Ha). If 
one now joins several segments together to form a long channel, then the radiation is asso- 
ciated with the junction points (top of one filament and bottom of the adjoining 'llainentl, 
and the coni|H>site channel radiates from its “kinks”. Tlie tortuous channel is one with 
many kinks, and therefore has many sources of radiation. As a result there is the potential 
for more structure (i.e. variability) in the signal radiated from tb»- tortuous channel than 
from the long straight channel which radiates only fron ,is two ends. Consequently, one 
would ex|x*ct that the spectrum of the tortuous chi' .nel will have its high frequency power 
increased in comparison with the equivalent straight channel, a conclusion which is sup- 
ported by the analysis and examines to be* presented below. 
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Within the limitation* inherent in the l-'raunhoffer appro ximation. either the “kink” 
pictM.e or the more conventional interpretation which ataociates radiation* with the fila- 
:nent center, yield equivalent re*ult*. Which picture one adopt* tlc|>end* on which of the 
two form* i* adopted for l(i') in Equation 4. For purpow* of the analyai* to follow, radi- 
ation will be a*aociated with the filament center; however, in the computer aimulation it 
wa* more convenient to aeparate the two e\ponential* in the *inc(x) function. 
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Mjllu’iiulically llic rjilulion uvit h> an obwncr ap|var» to he a »ct|iK*iicc of piiUe« 
emillcJ from each element with a temiHtral hbtory Ueivndent on the manner in which 
the puIh* pro|>a|(atek up the channel and on the luture (lenitth and orientatum) of the 
mdoidiial elementk of which the chaniH‘1 n c^impriH-d. AuKKiatinii thew puIm‘% with the 
filament center, the electric field emitted from the n*th clement a* seen b> an ohserser 
(for simplicity) on the siirfacv is obtained from I qiiation 


lutr.i') > J ^ - fU')a„sin('jkrfl^)iJ^*^'" 

»»} 


(<) 


w here 


^n»Pc„) 


A 

/ 




«»a) 


* 7 * ‘si j 

where I „ is the length ol the n-th element and in order to simplify the algebra, it has 
been avsumed that ‘VPcn •‘'••‘I field. I (r.»’). seen by the observer is the sum 

over all n of the I „(f.i'), and the spectrum of this radiation is defined to be: 


S(i>) 


v^<r(ra')r*(ra>» 


where the roinled brackets < > denote a statistical average. I'sing l-i|uatioii 5. one i^blaiiis 
the following form for the sivctrum: 
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|f(i/)|l ^ a„ sin('4kr}L„)ei*‘'^» 

1 * ^ 

^ an,sinC^kTjL,„)c j'“='m 

L w rj J iU « 


•Urn 


r M «n C4kTjL„)sin(VikTjLn)oi*‘'l^""^'" 1 

L irij J .Un 

•II m 


(7) 


In order to perform the required averages assume that: I ) The a„, 1^ and r„ are independ- 
ent random variables; 2) Tlie phase kefr^-rn,) is uniformly distributed over 2ir; and 3) 
a Pc where p^ is a constant, so that the a„ may be assumed to be identically distrib- 
uted. Then, the diagonal terms in the sum are dominant (e.g. Kodis, 1966) and one 
obtains: 

12 

|f(*»)| N<a„2><8in2('/,kTjL„)> (8) 

where N is tlie number of elements in the channel. In order to perform the average over 
the element lengths, L„, it will be assumed that the L„ are exponentially distributed with 
mean Lq. Based on observations of step|K*d leaders one would expect I.,, to be on the 
order of several lu's of meters (Uman, 1969; Schonland, 1956); however, little information 
exists conc'erning the actual statistics of the L^. Fortunately, the limiting behavior of the 
average at higli and low frequencies, which is the im|H>rtant factor for the discussion to 
follow, is independent of the choice of the distribution (Appendix A) and so the conclu- 
sions may be representative even if the assumed distribution is not. Also, the numerical 
results to be presented later support the conclusions reached here, but do not assume expo- 
nentially distributed L„. Performing the average over the L„ in equation 8, one obtains: 

<sin2('/2kf?Ln)> « .s (•>) 

I ‘♦•(KTyLQ r 


<E(r>)E*(r,i/)> a 


VpTT 
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I niployinit this roHiilt in l-iiuation N, one ohluin<( the lollowinK I'orin for the s|vetrum of 
rutliation from the torlmnu channel; 


S(»») 


N) 


I I +(krjly)* 




V4 


IftHI 


tIO) 


where 


S) 


»nPc 


|':N<a„2>|»4 


Notice that a* hitrh frequeiicies (kr)lQ»l ), S(i') is proportional to |f(i>)| and that at 
low fret|iK’ncies (krjl.o«l). •Sir') is prop-.riional to i»|f(i')|. Iliat is, in iVe hijih fre»|iieiK’y 
limit, the s|H'ctrum of radiated electric Held is proportional to the spectrum of the current 
pulsi*, whereas in the low friHiuency limit it is proportional to the spectrum of the current 
pulse multiplied hy frequency. I'he detlnition of low frei|uenc> must Iv treated carefull> 
here because for low enoui’h frequency the assumptions inherent in the I'raunhofler ap- 
proximation fail and the analysis is not correct. More precisely, the low freipiency behav- 
ior of S(i^) is to lx* re)!arded as the sha|x* of the spectrum in an intermediate frequency 
raiifse between (iou)!hly) krjLo * I and the low frequency bound on the l•■raunhoffer ap- 
proximation, roughly at kp^ * I. In this intermediate raniie S(t>) is proportional to i>|f(r')| 
' .ul at hi(!li frequencies Str*) is proportional to |f(i^)|. 

To appreciate the significance of these two regions, and to compare with data, the 
form of the current pulse is required. A pulse shape representative of currents measured 
in lightning return strokes is; 


l(t) = I, |e-“' -e-^^'l + Is le-"^' -e-*'| 


(II) 
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when* for typical firil return strokes a "2 x 10^, 0*8 x I0 *, 7« 10-' :ind 6 ■ 2 x lt>*, 
ami 1 1 ■ .U) kA and h " 2.5 kA. The shape of this cunent pulse is shown at the top in 
l-ilture 2 (ar.d its s|vctrum |f(i'>| is plotted at the hottoml. 11te first two exponentials in 
this expression represent the main current pulse in a form proposed by Bruce and (iolde 
(l')4l) with parameters sugfcested by IX*nnis and Piera* (IdM). The third term, h ex|><-7t), 
represent' intermediate current (Uman, 1^6*)); and the remaining exponential, h exp(> U). 
has been added to achieve continuity at t ■ 0 and doer^n'l, otherwise, significantly alter 
the current pulse. Details of the propagation of current during return strokes are still at 
issue (e.g.. Price and Pierce, I*)??). However, the preceding wave'orm reasonably represents 
a coni|H>site of reported current measuriMnents (liman, |dpd; ix*nnis and Pierce, ld(>4). 

The Courier transform of the pulse given in l-quation 1 1 is: 


f(i») 


(2iri')‘ 


0 - a 


11 *1 


•111 + 


7^ 


ii — ^ 
lo n^JT^in 


Tirri 


(12) 


il(i')i IS plotted in Figure 2 (bottom) for the parameters listed above. Notice in particular, 

**• > 

that in the limit of large fiequency (i> » ^I2n), the function Iftr*)! decreases as and 
that If (i')| IS constant in the low t're«tiK*ncy limit. 

Now substituting I'quation 12 for |f U')| into I quation 10, one can examine the liigli 
aiul '*low" frei|uency Ivhavior of the spectrum. .Siiuv the sp**ctruni S(i») of the radiated 
fields IS proportional to |l (»')| in the higli frequency him!, S(i') decreases as I /r‘ at high 
enough frequencies independent of toituousity. However, at intermediate frequencies - 
large enough that > 01 2n but small enough kijlg < I — the s|H*ctrum S(i') is propor- 
tional lo »'!'i(r')| and therefore decreases only as \/v. This slower rale of decrease at inter- 
mediate frequencies is a direct consequence of lortuousily. In contrast, the s|vclrum of a 


II 
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»in|{lo straighl channel (no tortuoustily ) not exhibit an inlermeiliate region with l/i^ 
liccrcaik* but rattier reaches a |vak and then rolls olT as |I'(»'H which is to say as l/v (Le 
Vine and Meneithini, l‘>7(«), I'xistence ol the intermediate reition lor the tortuous channel 
is consistent with available intormation on litihlnin^ i Murn strokes, l or example, letting 
I j) ■ .to meters, which is typical of the length of steps during the stepped leader (Uiiian, 
l‘>b‘>), and letting i} • ^ which is reasonable for return stroke velocities (Unian, l‘)(»‘M, 
one obtaii'.s krjl.o * I at ■ 0.5 MM/, but with ^ ■ H x 10*, If (I'M has already begun to 
decrease as l/i'- at these frequencies. 

riiis section has demonstrated how randomness of the channel can introduce energy 
at intermediate frequencies and thereby slow the rate of decrease of the spectrum, l x* 
aniples and a comparison with experiment will In* presented in the following sv'Ction where 
it will be shown that this effect is consistent with data. 
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In Ihii kccliun the ellecU of (ortuousify will be illuHlraleil with s>|>eelra ealeulaleil Iroiii 
Himiilaleii liithlning channeN. A siimilalion lor the purinm* of %tiulying the elfeet<k of loitu- 
ousity has been ilevelo|H*d by this author and a colleague (R Meneghinil using a piecewise lin- 
ear nuHlel for the channel and employing l-\|uations 3 and 4 for the fields radiated fri>in each 
element (I e Vine and Meneghini, 1478b). In this simulation the channel is generated by plac- 
ing end-ti>-end linear segments chosen by the machine in accordance with assigned statistics, 
file prtK'ediire is suggestive of real cloud-to-ground lightning in which the channel is lormed 
in a series of discrete ste|>s called the step|>ed leader. In fact, the simulation has been written 
so that radiation from the stepfvd leader may also be studied. In inmt of the work done 
to dale, the channel elements were chosen by assuming that the change in the (cailesian) 
ciHirdinales rvH|uired to advance from one end point to the next of a given element were 
normally distributed. Typically, the x- and y-ciH>rd inales were idenlieally distributed with 
zero .nean and the /.-coordinate had non-zero mean and a standard deviation on the order 
of JO meters. Obviously, one ean amtrol the general direction of the channel through 
choice of the mean value ol the coordinate change, and can control the variations about 
this direction bv means of the assigned standard deviations. Typical element lengths in 
the channels usi’d here were 40 meters, but have ranged from about ,)0 meler> to .100 
meters, l-'igure .) shows a channel generated by the computer, in this case without branche 
(Hie simulation is three dimensional so the projection of the channel on orthogonal 
planes is shown.) Once the channel has been generated, the computer calculates both the 
temporal history of the electromagnetic fields produced by a current pulse propagating 


jIoiI); lhi!t k'lunnol uiul jImi IIu< nunilikK* of llu* I'ouniT lran%rorm («|H*clruinl of the radi- 
alctl waveform. Piis is done iisiiii! the M>liitions given in rt|ualions 3 and 4 and their lime 
domain e«|uivalent. and emplo>ing an assigned eurreni waveform such as given in I >|ualion 
1 1. I1ie maior eominitaliim.il proolem is keeping Iraek of the proivr phaw of the enreni 
in eaeh of the many elements whieh may be radiating at a given instant. (See I e Vine 
and Meneghim, P>78b for esamples and a eomparison in the lime domain of simulated 
waveforms with data.) 

I1ie inUnenee of lorinonsily on the sivelriim of the rauialed fields is illiisiraled in 
figure 4. for eomparison sake, the envelo|v of the speelrum. Slid, is plotted tor both a 
long straight ehannel (no tortiiousity) and an et|iiivalent lorluoiis channel. Ilie curreiU 
waveform given in I ijuation I I was iisenl with a velocity of pr«ipagation of v ■ c/1 in both 
cases. I'he tortuous channel used in this calculation is shown m figure and was gene* 
rated from segments whose mean length was about 40 meters with variance of each carte- 
sian coordinate ol about -0 meters. In this esample, the observer is 50 km away and is 
located on the conducting surfice. The length of the vertical channel laboul 5 km) was 
chosem Si) that its spectral pea) was nearly identical with that of the tortuous channel. 

Piis was dime t<' facilitate comparisim of sliapes of the two spectra. (See I e Vine and 
Meneghini ( l‘>7(>) for a discriptiim of the effects i)f other parameters such as channel 
length and orientatit)n. and veliicity of priipagation on the spectrum of the straight channel. | 

Doth spectra shown in figure 4 |vak near 10 kP/. but the spectrum associated with 
the straight channel falls off at high fretiuencies as • /r-' whereas the s|vctrum associated 
with the tortuiuis channel decreases more slowly, first roughly as I 'i». and then near abi)ut 
4 Mil/ Ivgins to decrease .is I i'^. Pie high frei|uency asymptote in Imth cases is deter- 
mined by the current waveform (as discussed above) and therefore should be I however. 
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the variability introduced by tortuousity intrtxluces additional energy at intermediate fre- 
qucncicf blowing the rate of decrcaw of the tpectrum in this l'rci|uency reginu'. These 
spectra have the sha(H* deduced from the analysis of the preceding section. 

Ilie spectrum associated with the tortuous channel is consistent with measured spec* 
tra. fliis is illustrated in h'igure S where representative data (Kimpiira. l^bS) have been 
plotted with the calculated s|H‘ctra shown in l-igure 4. (Tlie data points were obtained 
from h'igure I of Kimpara, 1965, and normalised to unit bandwidth and an observer SO 
km from the channel.) As can be seen the spectrum associated with the tortuous channel 
(solid line) has a sha|H* consistent with the data n ; intermediate freitucncies; in c'ontrast. 
the $|K‘ctrum of the equivalent straight channel fulls off too rap'dly to agree well with 
the data, hventually, at a frequency determined by the mean length of the elements 
in the channel, the s|H'ctrum computed from the tortuous channel also begins to decrease 
as l/i/‘. llie transition from Mv to l/r»* decay is not inconsistent with data at high fre- 
quencies (e.g. hercc, 1977) and in fact Oh (Oh, Idod) interprets the data to have such 
a transition, although in the vicinity of 100 MM/. If confirmed by measurements, such 
a transition would have important implications. I'irst, establishing the frequency 
regime in which the transition occurs would provide a measure of the scale of chan- 
nel tortuousity (i.e. the mean element length) and by inference of the step length in 
the stepiH'd leader. Secxmdly, establishing the rate of decrease of the spectrum at frequen- 
cies beyond the transition would yield insight into the shape of the current pulse because 
in this regime the spc'ctrum is dominated by the pulse shape. Certainly, if a transition to 

l/r-2 

were to be found in the data, it would support the transmission line model and the cur- 
rent waveform in most common use (I'quation 10). 
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Sl'MMARY 


llitf liciuTil slw|V •»! Iho »|H*ctriiin ohiaiiicil hy iiuliulinic oHocI* i>l l4*rluouMly in 
kivpmit with ilala: It protlicU a |vak in the corrcii IrotiiioiK-y roKimo lollowal b> a ilo- 
iTcaso a» llii* inverse of fre»jm**'vy, which ihtmsIs Io a fret|uency depeiulenl on ihe lenulli 
i»f Ihe eleinenis in Ihe channel. Iliis s|H*clriim is in heller agri*einenl wilh ilala ihan Ihe 
s|vclriim pretliclcd hy Ihe name imHiel wilhoiil lorliioiisily (i.e, Ihe lonp siraiithi chaiinell. 
which decays al a rale proporlnnial lo I >^. 

Iliesi’ conclusions were derived frtvin a piecewisi*-linear ‘‘Iransinission line" iiuhIcI lor 
Ihe channel holh usiiift siinphlyinit assiimplions lo arrive al ch»vd form expressions for 
Ihe s|Hviriiin and fr»>ni numerical calciilalions on simiilaled channels. Hie anal>sis sii^esls 
an evenlual Iraiisilion h' I /r* deae.ise even for Ihe lorluous channel, which sui;t!esls a 
|H>lenlial means for leslintt Ihe model and itaining Insi^hl inlo relurn siroke ciirrenl and 
channel paramelers. 
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m.l KI CANIONS 

i'lieuri' I. Siimiiury ol im'>ikur«.'iiK*nl> I'roiii lighliimii. nicM.' iIjU Kjvc (hh'I iuh- 

nuli/i'il lo d tliMdiKV ol' 50 km diiJ unit hanilvkullh. Hu* <|H*clrum i» inoporliondl to 
irritiiK'iicy r* m llu' rriHiiKMU'y raiiitc I'rom j few lO't of kiloherl/ lo jboul 100 
nu‘)!Jherl/. 

|j|ture Hie eurreiil piilw uwd in llte edleuljlion« diiii IN <i|H\'triini li.e. nugniluile of 
lU I ourh'r irjiiklorniK llu* eurrenl pulw i< repreiu‘ntdli>e of current in fir<»l return 
kiroket. In Ihi* example, o ■ 2A ■ 4 x UH kec "* : • K x tO^ »i*c*"* ; > • 10-' ki*c"* 

and Ij) ■ 40 k.\ and l| • -.5 kA. 

I'lliure 4. I xampte of a lorluouii eliannel prodiuvd h> the computer ximulatum. 

l-igure 4. Ilie \|H*clrum calculated from a klrai|chl channel (about 5 km long) and a tor- 
tuous clianr i (the one slumn m I'lgure 4). Only the envelopes of the calculated 
s|vctra are pUdted. Ilie current pube shown in I’lgure ' was employed assuming a 
vehvity of propagation eipial to c/2. Flu* obsi*rver is 50 km from the channel. 

riguie 5. The caleulated spectra (from I'lgure 4) with represi-nlative data points su|vr- 
imposed. Pie data were taken from Kimpara (Kimpara, I^(v5) normah/ed lo 50 km 
and unit bandwidth. 


:o 



Figure I: Summary of spectral measurements from lightning. These data have been norinalued to a dKUnce of 50 km and 
unit bandwidth. The spectrum is proportional to ffrequencyr' in the frequency range from a few lOs of luloherla to about 
100 megahert/. 



FREQUENCY(kHi) 


I'il.uro 2 : Hio curront puK* uscti in tho calculations anil its spcciriini (i.c. inattniliulc of 
Its l ouricr transform). Hic current pulse is repress'iitative of current in first irfu; a strokes. 
In this evainple.tt = • 4 x l()^ sec"* ,/J * « x 10^ sec"' . > « 10-' sec"' a» ' In « 40 k \ 











(qp U| IAI/S-A) IAinU133dS 
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Figure 5. The calculated spectra (from Figure 4i with representative data points superimposed. The data were taken from 
Kimpara (Kimpara, 1965) normalized to 50 km and unit bandwidth. 



APPI NDIX A 


In this appcnilix the average <sin^ (ax)> will be performed in the limit of large and 
small a. It will be shown that the limits are independent of the probability density tunc* 
tion, Kx). lor a large class of hx). 


A. The large a limit: 


<sin2fax)> 


f 


P(x)sin2(ax)dx 


0 

which is absolutely convergent because |P(x)|dx ■ I. Now letting sin^(ax) ■ VA\ - 
cos2ax), one obtains: 


<sin2(ax» 


- 


f 


Ifx) eo8(2ax)dx 


But the integral on the right is zero in the limit o -• oo for piecewise continuous P(x) (c.g. 
see the Riemann-I.esbet|ue theorem, or Churchill, 1969). 


B. The small a limit: 

I'rom the theory for power series *)iie can write: 

A v^ 

sin2(y) * y2 + (1 - 2 cos^y)-^ 

where the last term on the right is the remainder after 4 terms of the Maclaurin series for 
sin^(y) and y is some particular y chosen to make both sides equal (c.g. Kaplan, 1952). 
With this result one obtains: 


A-1 


'f' 

<sin-(ax» ■ a2 j x’h x)dx f ( I -icos^y) h xklx 

0 0 

■ 0 “<x 2 > + -J I -2co»2y) <x^> 


Now noiinj! that |l-2cos“y| < I. and assuminic that the second and fourth moments of x 
exist, one has for sufficiently small a: 


sin-tttx) 2 a-<X“> 


Summari/ini!, <sin“(ox)> approaches Vi for very large a and is proportional to a- for 
Milficicntly small a. A s|H.’clfic example has been calcula.ed in the text emplo>ing expo- 
nentially distributed x with mean for which one obtains 

^ . S ^ 1 

<sin*(ax)> * 

“ 1 +(2oL^))’ 

I his is clearly equal to 1/i for very large o and proport ic/nal to a* for small a. 
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